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Abstract. We have undertaken high-resolution follow-up of a sample of high velocity H i clouds apparently associated with 
M31. Our sample was chosen from the population of high-velocity clouds (HVCs) detected out to 50 kpc projected radius of the 
Andromeda Galaxy by Thilker et al. (2004) with the Green Bank Telescope. Nine pointings were observed with the Westerbork 
Synthesis Radio Telescope to determine the physical parameters of these objects and to find clues to their origin. One additional 
pointing was directed at a similar object near M33. 

At 2' resolution we detect 16 individual HVCs around M31 and 1 HVC near M33 with typical Hi masses of a few times 
10^ Mq and sizes of the order of 1 kpc. Estimates of the dynamical and virial masses of some of the HVCs indicate that they 
are likely gravitationally dominated by additional mass components such as dark matter or ionised gas. Twelve of the clouds 
are concentrated in an area of only 1° x 1° at a projected separation of less than 15 kpc from the disk of M31. This HVC 
complex has a rather complicated morphological and kinematical structure and partly overlaps with the giant stellar stream of 
M31, suggesting a tidal origin. Another detected feature is in close proximity, in both position and velocity, with NGC 205, 
perhaps also indicative of tidal processes. Other HVCs in our survey are isolated and might represent primordial, dark-matter 
dominated clouds. 

Key words. ISM: clouds - Galaxies: Local Group - Galaxies: individual: M31, M33 - Galaxies: ISM - Galaxies: evolution - 
Cosmology: dark matter 


1. Introduction 

High-velocity clouds (HVCs) were discovered by Muller et al. 
(1963) in the 21-cm line emission of neutral atomic hydrogen. 
They are identified by their peculiar radial velocities which are 
incompatible with participation in the normal rotation of our 
Galaxy. HVCs can be found all across the sky in a variety 
of different sizes and shapes (Wakker & van Woerden 1997). 
Some HVCs form large and coherent complexes spanning tens 
of degrees in size. Others are comparatively small and isolated. 
The origin of HVCs is still under debate. In the decades since 
their discovery many different mechanism have been proposed 
to explain the observed population of HVCs (see Wakker & van 
Woerden 1997 for a detailed review). One of the most promi¬ 
nent hypotheses puts forward the Galactic fountain (Shapiro 
& Field 1976) as a means of forming HVCs. In this scenario, 
hot gas ejected by the Galactic fountain would cool down and 
recombine high above the Galactic plane. The condensed gas 
clouds, falling back ballistically towards the Galactic plane, 
would then appear as neutral high-velocity clouds. Other sce¬ 
narios suggest that the gas of HVCs was torn out of dwarf 
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galaxies during close encounters with the Milky Way by tidal 
forces or by ram-pressure interaction. At least two promi¬ 
nent structures classified as HVCs, the Magellanic Stream 
(Mathewson et al. 1974) and the Leading Arm (Gardiner & 
Noguchi 1996; Putman et al. 1998), were formed in this way. 

Another hypothesis that has been discussed extensively in 
recent years postulates that HVCs are primordial gas clouds 
left over from the formation of the galaxies in the Local Group. 
This scenario was first suggested by Oort (1966) and recently 
revived in the framework of cold dark matter (COM) cosmol¬ 
ogy by Blitz et al. (1999), who found the kinematic proper¬ 
ties of most HVCs to be consistent with a distribution through¬ 
out the entire Local Group. They proposed that HVCs might 
represent the gaseous counterparts of the so-called “missing” 
dark matter satellites which were predicted by CDM structure 
formation scenarios (Klypin et al. 1999; Moore et al. 1999; 
Kravtsov et al. 2004). On the basis of the Leiden/Dwingeloo 
Survey (LDS) of Galactic neutral hydrogen (Hartmann & 
Burton 1997) Braun & Burton (1999) identified a sub-sample 
of 66 compact high-velocity clouds (CHVCs) with angular 
sizes of less than 2°. The kinematic parameters of the CHVCs 
suggested an unvirialized population within the Local Group 
potential with distances that might extend to the turn-around 
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Right Ascension (J2000) 


Fig. 1. Overview of the 9 fields (hlack circles, labelled with cap¬ 
ital letters) observed around M31 with the WSRT. The grayscale 
image displays the Hi emission of M31 observed by Thilker et al. 
(2004) with the GET. The contours show our Hi data of the 16 high- 
velocity clouds detected with the WSRT (labelled with numbers) and 
of NGC 205. In fields H and I no reliable signal of high-velocity Hi 
emission was found at the resolution and sensitivity of our WSRT data. 



Right Ascension (J2000) 


Fig. 2. The map compares the position and heliocentric velocitiy of 
the giant stellar stream of M31 with the nearby high-velocity clouds. 
The squares indicate the positions of the radial velocity measurements 
along the stream (outlined by the dotted lines) by Ferguson et al. 
(2004). For comparison we have inserted the first moment maps of 
the HVCs from our data. The GET data of the Hi disk of M31 from 
Thilker et al. (2004) is shown for orientation in black contours in the 
upper right corner of the map. The numbers denote the heliocentric 
radial velocities of individual features. 


radius of about 1 Mpc, representing the primordial dark-matter 
mini-haloes which have not yet been accreted by one of the 
large galaxies. 

Recent observations have raised doubt on an extended 
Local Group scenario. If the observed CHVC population had a 
typical distance as large as 1 Mpc, implying a typical H i mass 
of 10^ Mq, then we would expect to detect a similar population 
of Hi clouds in other galaxy groups. However, all attempts to 
detect such an intra-group population in nearby galaxy groups, 
employing H i mass limits of about 10^ Mq, have failed (Zwaan 
2001; Braun &. Burton 2001; Pisano et al. 2004). In addition, 
Ha observations towards a number of CHVCs (Tufte et al. 
2002; Putman et al. 2003) have been interpreted to imply dis¬ 
tances of the order of only 10 kpc from the Galaxy, where the 
Galactic radiation field is strong enough to account for the de¬ 
tected Ha intensities. All these observations suggest that the 
CHVCs are distributed throughout the circum-Galactic envi¬ 
ronment, more closely associated with the Milky Way rather 
than being concentrated near the Local Group turn-around ra¬ 
dius. 

The closest large spiral galaxy comparable to the Milky 
Way is the Andromeda galaxy (M31), making it an ideal tar¬ 
get to verify the circum-galactic nature of high-velocity clouds 
observationally. Thilker et al. (2004) used the Green Bank tele¬ 
scope (GBT) to observe a field of 7° x 7° around M31 in the 
21-cm line of neutral hydrogen. They discovered a population 
of Hi clouds (with Hi masses between 10^ and 10® M©, in the 


case of Davies’ Cloud about 10^ Mq) which are likely to be 
the counterparts of the high-velocity clouds around our Milky 
Way. A similar GBT survey of a 5° x 5° field around M33 by 
Thilker et al. (in prep.) covered another anomalous-velocity H i 
cloud close to the disk of M33. This cloud was already de¬ 
scribed by Huchtmeier (1978) and has properties similar to the 
HVCs found in the direction of M31. Many of the clouds in the 
direction of M31 and M33 are rather compact and barely re¬ 
solved with the GBT. Moreover, complex line profiles in some 
directions indicate the presence of compact substructure which 
is unresolved in the GBT data. Therefore, we re-observed some 
of the HVCs near M31 and M33 in Hi, using the Westerbork 
Synthesis Radio Telescope (WSRT). The aim of our interfero¬ 
metric observations was to study the HVCs and their compact 
substructure with higher angular resolution and to determine 
their physical parameters. 

In this paper we present the results of these WSRT obser¬ 
vations and discuss the various possible origins of the detected 
HVCs in the neighbourhood of M31 andM33. After describing 
the sample selection in Sect. 2, we specify the technical param¬ 
eters of the observations and the data reduction procedure in 
Sect. 3. In Sect. 4 the results of our observations are presented. 
A discussion of our results and their implications can be found 
in Sect. 5. Finally, Sect. 6 gives a summary of our results and 
conclusions. 
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Table 1. Physical parameters of the high-velocity clouds found near M31 and M33. The columns give the cloud name, the right ascension a 
and declination S, the peak column density Nhi at 2' resolution, the heliocentric radial velocity Vhd, the average line width (FWHM) of the Hi 
lines Av, the cloud diameter D (FWHM), the H i mass Mhi, the H i mass in the GET data of Thilker et al. (2004) the dynamical mass 

Mdyn, the virial mass Myir, the average H i volume density mhi, the value of the velocity gradient 6v, the position angle of the velocity gradient ip 
(measured in the mathematically negative sense with 0° being north), and the ellipticity e. Distances of 780 kpc for M31 and 800 kpc for M33 
were assumed for the clouds to calculate their sizes, masses, and densities. In the last two rows we have calculated the mean value {x) and the 
standard deviation of some parameters of M31 HVCs 1-16. 
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M31 HVC 1 

0''42"'53* 

38°52' 

11.1 

-503 

17.3 

0.72 

5.2 

4.8 

>86 

230 

10.9 

7.7 

349 

0.24 

M31 HVC 2 

0'’41”03" 
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9.5 

-511 

29.6 

0.75 

5.0 


>44 
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9.3 

9.7 
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0.12 

M31 HVC 3 

OMO^IH 

39°21' 

3.4 

-509 

18.5 

0.93 

2.2 
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2.1 

8.4 
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0.28 

M31 HVC 4 

0H3“07* 

39°20' 

4.2 

-518 

21.0 

0.56 

1.3 
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4.1 

98 

0.22 

M31 HVCS 
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0.48 
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0''42“3H 

39°41' 
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-517 

17.3 

0.53 
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3.2 

4.1 
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0.42 
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39°26' 

3.7 

-464 

20.3 

1.27 

3.5 




1.3 

10.7 

280 

0.40 

M31 HVC 8 

0'’41“47'‘ 

39°28' 

3.0 

-476 

32.1 

0.92 

2.1 


>730 


2.1 

15.5 

291 

0.36 

M31 HVC 9 

0'’43“13* 

39°39' 

1.8 
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19.1 
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0.8 



260 

1.9 

12.7 
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0.51 

M31 HVC 10 

0'’41“03" 

39°38' 
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30.2 

1.20 

5.0 
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2.3 

14.7 
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0.38 

M31 HVC 11 

0 H 3 “ 3 y 

39027/ 

0.9 

-454 

30.2 

0.43 

0.2 
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2.0 



0.26 

M31 HVC 12 

01'43“49* 

39°33' 

1.3 

-390 

30.8 

1.5 

1.5 




0.3 




M31 HVC 13 

0'’39“35* 

41°15' 

1.3 
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3.5 
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0.65 
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3.0 
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3.3 

4.5 
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0.37 
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25.9 
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3.9 
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0.43 
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0.57 

M33 HVC 1 
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5.2 
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2.6 



0.43 
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2. Sample selection 

Based on the GBT observations of Thilker et al. (2004) and 
Thilker et al. (in prep.) we selected 9 fields around M31 and 
1 field near M33 for follow-up observations of HVCs with the 
WSRT. The 9 fields around M31 were centred on the bright¬ 
est and most conspicuous clouds and cover about half of the 
population of HVCs discovered by Thilker et al. (2004). The 
targets cover the different morphological types found in the 
GBT data and include HVCs with both large and small an¬ 
gular separations from the disk of M31. The most prominent 
cloud near M31, Davies’ Cloud, was already studied in great 
detail (Davies 1975; de Heij et al. 2002) so we did not include 
it in our survey. An overview of the 9 fields around M31 can be 
found in Fig. 1. 

The field near M33 is located about 1° south of the centre 
of the galaxy. It is centred on a faint H i cloud outside the disk 
of M3 3 which is prominent in the GBT data. The cloud is not 
isolated but appears to be connected with the Hi disk of M33 
by a faint gas bridge. 

3. Observations & data reduction 

The 10 selected fields were observed in the 21-cm line emis¬ 
sion of Hi, using the WSRT in the Maxi-Short configuration. 


This configuration provides particularly good sampling of the 
shortest East-West baselines (36, 54,72 and 90 m). A complete 
12-hour track for each field was preceded and followed by ob¬ 
servations of the external calibrator sources CTD93 and 3C147. 
The correlator provided 1024 spectral channels for each of the 
two linear polarisations. To improve the signal-to-noise ratio, 
the 1024 channels were equally divided among two indepen¬ 
dent IF bands tuned to the same central frequency, resulting in 
512 channels for the final spectrum. By later averaging these 
two IF bands we are able to decrease the digitization noise by 
a factor V2. As a result, the noise level in the final images is 
about 5% lower in comparison with the normal single-IF mode. 
With a bandwidth of 5 MHz the corresponding velocity resolu¬ 
tion is about 2 km s^*. 

The data reduction was done with the AIPS software pack¬ 
age. After flagging the data affected by radio frequency inter¬ 
ference or by shadowing of the telescopes we carried out the 
standard bandpass, gain, and flux calibration, using the two ex¬ 
ternal calibrators. To further improve the gain calibration, we 
self-calibrated on the continuum sources in each field in an it¬ 
erative procedure. This provided us with a clean component 
model of the continuum emission which was then subtracted 
from the visibility data. 

As the emission of most HVCs in our survey was quite faint 
and extended, we applied a Gaussian uv taper, declining to 30% 
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Fig.3. (a) Hi volume densities «hi and (b) heliocentric radial velocities Vhei of M31 HVCs 1-12. M31 HVCs 1-6 (•) have similar radial 
velocities, and they are much more condensed with typically higher densities. M31 HVCs 7-12 (o), in contrast, are more diffuse with lower 
densities and rather different radial velocities. In (c) we have plotted the value of the velocity gradient across some of the clouds (labelled with 
italic numbers) in kms“* arcmin^' (labelled from 1 to 4) against the direction of the gradient in degrees. It turns out that the gradients across 
M31 HVCs 1-6 (solid arrows) have arbitrary orientations whereas the typically larger gradients across M31 HVCs 7-10 (dashed arrows) are 
all oriented in about the same direction. Note that for M31 HVCs 5, 11, and 12 we could not obtain a meaningful fit for the velocity gradient. 


amplitude at 1.25 kd radius, before producing the final, cleaned 
images. The resulting synthesised beam has a FWHM of about 
2', corresponding to a spatial resolution of about 450 pc at the 
distance of M31. For the brightest and most compact clouds we 
also produced images with a higher angular resolution of 30" 
FWHM (using a uv taper radius of 5 kd). The velocity reso¬ 
lution was also smoothed to 12 km s * in all spectra to further 
increase the signal-to-noise ratio, given that the typical detected 
linewidth was about 25 kms“^ For the 2' resolution maps, at 
12 km s^*, we obtain an RMS sensitivity of about 1 mJy/beam 
towards the centre of the field, corresponding to a brightness 
temperature RMS of 40 mK. This implies a Icr H i column den¬ 
sity sensitivity of 9 ■ 10'^ cm“^ for emission filling the 2' beam. 


4. Results 

In two of the nine fields around M31 no high-velocity gas was 
detected. In both cases the emission of the high-velocity gas is 
probably too faint and/or diffuse to be detected in our WSRT 
data given the finite column density sensitivity noted above (9 ■ 
10^^ cm“^ over 12 kms^*). In all other fields we clearly detect 
the HVCs at more than 5cr significance. The individual clouds 
have been named Mxy HVC z, where xy is the Messier number 
of the host galaxy and z is the catalogue number assigned to 
each HVC. 


The observational parameters of all clouds are summarised 
in Table 1. Under the assumption of a constant distance of 
780 kpc (Stanek & Garnavich 1998) we obtain a mean linear 
size of the 16 HVCs around M31 of about 1 kpc FWHM and 
a mean Hi mass of 3 ■ 10^ Mq. For the HVC near M33 we as¬ 
sumed a distance of 800 kpc (Lee et al. 2002) to determine the 
physical parameters. 

The HI column densities were derived from the zeroth mo¬ 
ment of the spectra under the assumption that the optical depth 
of the gas is negligible. From the integrated column density 
we can directly calculate the Hi mass of each cloud via 
Mhi = where imh is the mass of a hydrogen 

atom, d is the distance of the cloud, and tp denotes the angu¬ 
lar size of a resolution element of the map. The observed H i 
masses are in the range of a few times 10“* Mq to 6 ■ 10^ Mq. 
For the isolated clouds in the direction of M31 we have also 
listed in Table 1 the H i masses derived from the GET data of 
Thilker et al. (2004). These values are in good agreement with 
our WSRT observations. Only for the most extended and dif¬ 
fuse cloud, M31 HVC 13, we detect only about 40% of the total 
flux. 

The diameters of the HVCs were determined by fitting a 
Gaussian to the radial column density distribution. By divid¬ 
ing the H i mass by the volume defined by the FWHM of this 
Gaussian, we can estimate a mean H i volume density hhi for 
each cloud under the additional assumption of spherical sym- 
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Fig. 4. M31 HVC 1. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm 5 • lO'* cm and from 
1 • 10'^ cm“^ in steps of 1 • 10*^ cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“'. The 
diagrams on the left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of 
the gas along the cut marked hy the arrow in the column density map. 


metry. The derived densities are typically of the order of a 
few times 10'^ cm~^. These volume densities are between 
two and three orders of magnitude higher than the densities 
of 10“^ ... 10 cm“^ expected for an extended galactic corona 
in which the HVCs might be embedded (Sembach et al. 2003; 
Rasmussen et al. 2003). 

The mean radial velocities of the HVCs were calculated 
from the first moment of the spectra. Most values are more neg¬ 
ative than the M31 systemic velocity of Vhei ~ -300 kms“*. 
However, this is a selection effect, because most of the HVCs 
detected around M31 belong to a filamentary complex of 
clouds with similar radial velocities near the southern edge of 


the galaxy. The HVC near M33 has a slightly less negative 
radial velocity with respect to the M33 systemic velocity of 
Vhei ~ -180 kms“*. By fitting a Gaussian to the spectral lines 
we extracted the line widths with typical values in the range of 
about 20 .. .30 kms^* FWHM at 2' resolution. These values 
correspond to an upper limit for the kinetic temperature of the 
gas of the order of 10"^ K. Similar line widths are observed for 
the Galactic HVCs (e.g. Wakker & van Woerden 1997, Braun 
& Burton 1999, Westmeier et al. 2005b), and they are typi¬ 
cal for the warm neutral medium (Braun 1997). Thilker et al. 
(2004) already noted that the large line widths of the HVCs 
discovered in the GBT survey of M31 require a large fraction 
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Right Ascension (J2000) 

Fig. 5. Map of the H i column densities (top) and the heliocentric ra¬ 
dial velocities (bottom) of M31 HVCs 2-12. The contours in the upper 
map are drawn at 1 • lO'* cm“^, 5 • lO'* cm“^ and from 1 • lO'^ cm“^ 
in steps of 1 • 10*^ cm“^. The contours in the lower map are separated 
by 5 km s“'. 


of undetected mass to allow for gravitational stability of the 
clouds. We address this problem in more detail in the discus¬ 
sion in Sect. 5.1. 

We also determined the radial velocity gradients across 
some of the HVCs by fitting a plane P(x,y) = cx + dy + f 
to the first moment maps Vh^\{x,y). The chi-square minimisa¬ 
tion provided us with values for the parameters c, d and / from 
which we calculated the absolute value 6v - and the 

direction (p — 90° - arctan(t//c) of the velocity gradient. In 
this definition the position angle, p, is measured in the usual 
(east of north) sense. The velocity gradients are listed in Table 1 
for those clouds where a meaningful fit could be obtained. The 
method, of course, only provides us with a linear approxima¬ 
tion of an average global velocity gradient across each cloud, 
whereas the observed gradients are not expected to be perfectly 
linear and, on the local scale, will deviate from the average 
value and direction determined from the planar fits. 

The ellipticities of the HVCs were obtained by fitting an el¬ 
lipse to the column density distribution. The parameters of the 
ellipse were derived from a weighted second moment analysis, 
following the procedure described by Banks et al. (1995). The 
ellipse centre was fixed by the first moment of the column den¬ 
sity distribution. From the calculation of the second moments 
we then derived the major and minor axis, a and b, of the el¬ 


lipse. The ellipticities given in Table 1 are defined by e = 1 - 
Only M31 HVC 12 was so diffuse and extended that a mean¬ 
ingful fit could not be obtained. 

4.1. M31 HVC 1 

M31 HVC 1 is located about half a degree south of 
M31 HVCs 2-12. This proximity together with similar ra¬ 
dial velocities suggests a connection with M31 HVCs 2-12. 
M31 HVC 1 reveals the highest peak column density detected 
among all clouds in our sample with Nhi = 1.1 ■ 10^*^ cm“^ 
at 2' resolution. The Hi mass is Mui = 5 ■ 10^Mq. Maps 
and profiles of the cloud are presented in Fig. 4. At first 
glance, M31 HVC 1 has a spherically-symmetric appearance. 
Furthermore, it shows a systematic radial velocity gradient, 
measuring about 15 km s ' in an approximately north-south 
direction. The high brightness of M31 HVC 1 allowed us to 
additionally generate a map with 30" angular resolution which 
is shown in Fig. 15. The high-resolution column density map 
resolves M31 HVC 1 into a fairly compact core with a peak 
column density of more than 2 ■ 10^" cm“^ which is embedded 
in a more diffuse envelope with a conspicuous extension in the 
northern direction. We also produced a data cube with higher 
velocity resolution which shows that the spectral lines in the 
direction of the core have significantly smaller line widths of 
only about 6 km s ', indicating the presence of cold gas in the 
core. 

4.2. M31 HVCs 2-12 

M31 HVCs 2-12 are probably the most remarkable feature de¬ 
tected in our survey. They are crowded in an area in the sky of 
about 1° in size, and they are located in projection close to the 
edge of the H i disk of M31. An overview of the whole field 
is shown in Fig. 5 while maps and profiles of some individ¬ 
ual clouds are presented in Fig. 6, 7, 8, and 9. Maps with 30" 
resolution of some of the HVCs are shown in Fig. 15 and 16. 
The individual clouds in this complex of HVCs appear to be ar¬ 
ranged in filaments intersecting each other at different angles. 
The radial velocities of the clouds cover a very large range of 
-520 kms^' < Vhei ^ -390 kms^', and it looks as if clouds 
along the apparent filaments have similar radial velocities. 

There also appears to be a morphological distinction be¬ 
tween some of the filaments. M31 HVCs 2-6 (and also HVC 1 
which is presumably associated with the HVC complex) have 
very similar radial velocities of Vhei ~ -520... - 500 kms^'. 
Furthermore, they all show high peak fluxes and mean volume 
densities as well as being fairly condensed. M31 HVCs 7- 
12, in contrast, have different radial velocities of Vhei ~ 
-480... -390 km s“', and they are more diffuse with typically 
lower peak fluxes and densities. The average H i volume densi¬ 
ties and the heliocentric radial velocities of the 12 HVCs are 
plotted in Fig. 3 (a) and (b). A distinction between the two 
different groups of HVCs in this area can also be made on 
the basis of the velocity gradients observed across the clouds. 
While the velocity gradients across M31 HVCs 1-6 have ar¬ 
bitrary orientations, the somewhat larger gradients observed 
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Fig. 6. M31 HVC 2. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm 5 • lO'* cm and from 
1 • 10'^ cm“^ in steps of 1 • 10*^ cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“'. The 
diagrams on the left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of 
the gas along the cut marked hy the arrow in the column density map. 


across M31 HVCs 7-12 appear to be oriented in the same di¬ 
rection (see Fig. 3 (c)). 

Two other clouds at the north-western edge of the map have 
been named M31 HVC A and B. Both are fairly diffuse with 
only moderate H i peak column densities and relatively broad 
spectral lines which distinguishes them from the HVCs dis¬ 
cussed before. The radial velocities along the two clouds seem 
to follow those observed for the H i disk of M31 which was also 
covered at the edge of our field of view and which can partly 
be seen in the upper-right corner of the column density map in 
Fig. 5. This suggests that M31 HVC A and B might represent 
gas of an outer disk component of M31. 


4.3. M31 HVC 13 

M31 HVC 13 is rather diffuse and extended. It has quite low 
column densities and appears to be fragmented into smaller 
clumps towards the eastern edge. Maps and profiles are shown 
in Fig. 10. M31 HVC 13 appears in projection against the 
HI disk of M31. But its radial velocity is different by about 
130 km s * from that of the disk gas observed in this direction, 
so that it is presumably located outside the disk plane of M31. 
It lies only about half a degree south of NGC 205, the satellite 
galaxy of M31. H i emission from NGC 205 was also detected 
at the edge of our field of view. Interestingly, the radial veloc- 
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Fig.7. M31 HVC 3. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm“^ and from 5 • 10** cm“^ in 
steps of 5 • 10** cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“*. The diagrams on the 
left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the 
cut marked by the arrow in the column density map. 

ities of the gas observed in NGC 205 and M31 HVC 13 cover ance. While M31 HVC 14 is elongated and slightly bent, 
essentially the same range. Despite the positional and velocity M31 HVCs 15 and 16 exhibit a pronounced head-tail struc- 
correspondence we do not detect a continuous H i bridge be- ture. A compact core appears to be embedded in an asymmet- 
tween these two objects at our sensitivity. ric envelope which forms an extended tail in one direction. 

Head-tail structures like these are also widely observed in the 
case of Galactic HVCs (Briins et al. 2000; Briins et al. 2001; 
4.4. M31 HVCs 14—16 Westmeier et al. 2005b), and they are suggestive of a distortion 

of the clouds by external forces like e.g. the ram pressure of an 
The remaining three HVCs observed near M31 are all posi- ambient medium, 
tionally isolated from each other and from the disk of M31. 

Maps and profiles are shown in Fig. 11, 12, and 13. All three 
clouds clearly deviate from a spherically-symmetric appear- 
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Fig. 8. M31 HVC 4. The upper right map displays the Hi column density. The contours are drawn from 1 • lO'^ cm ^ in steps of 3 • lO'* cm 
The lower right map shows the heliocentric radial velocity. The contours are separated by 5 kms“*. The diagrams on the left hand side show 
from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the cut marked by the 
arrow in the column density map. 


4.5. M33 HVC 1 


5. Discussion 


M33 HVC 1 is located about 1° south of the centre of M33. It 
is the most conspicuous HVC around M33 found in the GET 
survey by Thilker et al. (in prep.). In the GET data it is already 
obvious that the HVC is not isolated but connected with the H i 
disk of M33 by a faint gas bridge. In our WSRT data (Fig. 14) 
this bridge is very prominent with H i column densities of a few 
times lO'® cm“^. Eoth M33 HVC 1 and the gas bridge share a 
common gradient in radial velocity measuring about 20 km s“*. 


5.1. Dynamical mass and stability 

Three clouds, M31 HVCs 1, 2, and 6, have a spherical appear¬ 
ance combined with a pronounced gradient in radial velocity. 
Two more clouds, M31 HVCs 8 and 10, have a pronounced, 
linear velocity gradient, although they appear slightly ellipti¬ 
cal. If we assume these velocity gradients to be caused by a 
rotation of the clouds we can estimate their dynamical mass 


^dyn ~ 
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Fig. 9. M31 HVC 6. The upper right map displays the Hi column density. The contours are drawn from 1 • lO'^ cm ^ in steps of 3 • lO'* cm 
The lower right map shows the heliocentric radial velocity. The contours are separated by 5 kms“*. The diagrams on the left hand side show 
from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the cut marked by the 
arrow in the column density map. 


where R is the radius of the cloud, Vrot is the rotation ve¬ 
locity at the edge, and G denotes the gravitational constant. 
The velocity gradients 6v over the angular extent t? were es¬ 
timated from the one-dimensional cuts across the clouds and 
are summarised in Table 2. Because of the unknown incli¬ 
nation angle i of the clouds the observed velocity gradient 
6v only gives us a lower limit of the rotation velocity with 
Vrot = (5v/(2cos/). Nonetheless, the obtained lower limits for 
the dynamical masses of the five clouds are already by one or 
two orders of magnitude larger than the observed H i masses 
(see Table 1). If the assumption of rotation is correct we are 
seeing only a small fraction of the clouds’ masses in neutral 


atomic hydrogen. Additional mass components like ionised gas 
or dark matter are required to account for the observed velocity 
gradients. 

A similar result can be derived by applying the virial the¬ 
orem, assuming spherical symmetry and, for simplicity, a con¬ 
stant mass density. We can justify the assumption of virialisa- 
tion by making a simple estimate of the dynamical timescales 
of the clouds. If we consider the mean diameter of the clouds 
to be D ss 1 kpc and the mean FWHM of the spectral lines 
to be Av ss 25 km s“* we obtain an estimate for the dynam¬ 
ical timescale of Tdyn = DjAv ^ 4 • 10^ yr. This timescale 
is short compared with the expected orbital timescales of the 
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Fig. 10. M31 HVC 13. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm“^ and from 5 • 10** cm“^ in 
steps of 5 • 10** cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“*. The diagrams on the 
left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the 
cut marked by the arrow in the column density map. 


HVCs relative to M31, about 10® yr, over which changing tidal 
forces might be active. Thus, we can assume that the HVCs are 
at least close to internal dynamical equilibrium so that we can 
apply the virial theorem which, in our case, reads 


Here, Mvir is the virial mass of the cloud, R is the cloud radius, 
G denotes the gravitational constant, and Av is the FWHM of 
the HI lines. Solving Eq. 2 for the virial mass and inserting the 
parameters given in Table 1 yields virial masses of the order of 
a few times 10® Mq for the most condensed clouds in our sam¬ 


ple. The calculated virial masses lead to a typical mass fraction 
of neutral hydrogen of less than 1%. Although these are small 
values, we have to consider that the observed line widths are 
probably only upper limits for the velocity dispersion of the 
gas since internal kinematics or turbulence can result in an ad¬ 
ditional broadening of the spectral lines. 

It should also be stressed that the application of the virial 
theorem to H i clouds is open to justified criticism. Not only 
does the calculation of virial masses depend upon a number of 
assumptions with respect to the object symmetry, but the sim¬ 
ple form of Eq. 2 also demands a homogeneous population of 
particles with identical masses. It is doubtful that the velocity 
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Fig. 11. M31 HVC 14. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm ^ and from 5 • 10** cm ^ in 
steps of 5 • 10** cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“*. The diagrams on the 
left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the 
cut marked by the arrow in the column density map. 


dispersion of the H i component alone would provide a realistic 
estimate of the total mass of HVCs. Moreover, the calculation 
of dynamical masses according to Eq. 1 is in conflict with the 
virial theorem in the form of Eq. 2. The virial theorem requires 
random motions of the particles instead of a regular rotation of 
the cloud. Either the dynamical mass or the virial mass of a par¬ 
ticular cloud can be correct but not both. By calculating virial 
masses, we can only expect to assess the order of magnitude of 
the entire mass of a cloud which, in our case, confirms the ne¬ 
cessity of a large amount of additional mass such as ionised gas 
or dark matter to stabilise the HVCs observed in our survey. 


5.2. The origin of the HVCs 

One possible origin of the HVCs observed around M31 and 
M33 could be tidal stripping of gas during close encounters 
with dwarf galaxies. Thilker et al. (2004) already noticed that 
the complex of HVCs found south of the disk of M31 is partly 
overlapping with the giant stellar stream of M31 discovered 
by Ibata et al. (2001). A comparison between the optical map 
of the stellar stream by Eerguson et al. (2004) and our WSRT 
Hi data (Eig. 2) shows that the HVCs are grouped along the 
southern edge of the stream close to the disk of M31. This 
positional overlap suggests a connection between the stellar 
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Fig. 12. M31 HVC 15. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm“^ and from 5 • 10** cm“^ in 
steps of 5 • 10** cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“*. The diagrams on the 
left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the 
cut marked by the arrow in the column density map. 


stream and the H i gas. Ferguson et al. (2004) also present the 
results of radial velocity measurements of the stellar compo¬ 
nent at four positions along the stream. The heliocentric radial 
velocities range from -550 km s“* close to the disk of M31 up 
to -325 kms“* about 4° away from the disk. In the region of 
the HI gas, the radial velocities of the stellar component are in 
the range of about -480 km s^^ to -550 km s^^ which is partly 
overlapping with the velocities found for the HVCs in our sur¬ 
vey. But, as noticed in Sect. 4.2, the velocity structure of the 
HVC complex is very complicated. Several filaments of clouds 
with different radial velocities appear to intersect each other at 


various angles, making a comparison with the velocities of the 
stellar stream challenging. 

In addition to having a degree of organization into filamen¬ 
tary structures, it is striking that many of the H i clumps, par¬ 
ticularly HVCs 2-6 seen in Fig. 5, have large internal veloc¬ 
ity gradients which are typically not directed along the fila¬ 
ments. This circumstance suggests that even within this region 
of likely overall tidal origin, there may well be dynamical ef¬ 
fects related to localised dark matter concentrations. 

Another example of a cloud that is likely of tidal ori¬ 
gin is M31 HVC 13 which is located about half a degree 
south of NGC 205 and covers the same radial velocities as the 
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Fig. 13. M31 HVC 16. The upper right map displays the Hi column density. The contours are drawn at 1 • lO'* cm“^ and from 5 • 10** cm“^ in 
steps of 5 • 10** cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“*. The diagrams on the 
left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the 
cut marked by the arrow in the column density map. 


gas observed towards this satellite galaxy of M31. A possible 
tidal interaction between NGC 205 and M31 was already sug¬ 
gested by Zwicky (1959), based on the discovery of extended, 
faint arms emanating from both ends of NGC 205. Recently, 
McConnachie et al. (2004b) reported the discovery of an arc¬ 
like feature in the distribution of blue red giant branch stars 
north of NGC 205 which might be part of a stellar stream orig¬ 
inating from a tidal interaction with M31. This presumable tidal 
trail, however, is located on the opposite side of NGC 205 with 
respect to M31 HVC 13. Our current sensitivity has not per¬ 
mitted detection of a continuous Hi bridge between the two 
objects, so a physical association remains circumstantial. 


We have already noted that M33 HVC 1 is not isolated but 
appears to be connected with the H i disk of M3 3 by a faint gas 
bridge. This possible connection is suggestive of a tidal ori¬ 
gin of M3 3 HVC 1 although the potential progenitor remains 
unkown. The distribution of stellar sources in the INT WFC 
survey of M33 (McConnachie et al. 2004a) appears to show 
no significant enhancement toward M33 HVC 1 or its associ¬ 
ated HI bridge. This suggets the progenitor of M33 HVC 1 ei¬ 
ther: (1) is quite low in stellar mass relative to the Andromeda 
stream, (2) lies outside the INT survey field, or (3) has been 
completely integrated by M33. It is also possible that a pro¬ 
genitor never existed, that is, M33 HVC 1 never supported star 
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Fig. 14. M33 HVC 1. The upper right map displays the Hi column density. The contours are drawn at 1 • 10** cm ^ and from 5 • 10** cm ^ in 
steps of 5 • 10** cm“^. The lower right map shows the heliocentric radial velocity. The contours are separated by 5 km s“’. The diagrams on the 
left hand side show from top to bottom: the H i column density, the heliocentric radial velocity, and the velocity dispersion of the gas along the 
cut marked hy the arrow in the column density map. 


formation and remains largely intact as a gaseous, primordial 
dark matter halo. Thus, the region around M33 HVC 1 remains 
a promising target for deep optical surveys which could un¬ 
cover peculiarities in the stellar distribution associated with the 
high-velocity H i gas. 

The remaining three HVCs studied near M31 are all iso¬ 
lated in position-velocity space from both the stellar stream 
and any known dwarf companion of M31, making a tidal origin 
questionable. Another possible origin for the HVCs arises from 
current CDM structure formation scenarios, predicting a hierar¬ 
chical formation of dark-matter haloes. Large galaxies like the 
Milky Way or M31 are believed to have formed by accreting 


smaller dark-matter haloes and thereby growing to their cur¬ 
rent mass and size. Numerical simulations predict a large num¬ 
ber of so-called dark-matter mini-haloes throughout the Local 
Group to be left over from this hierarchical process of galaxy 
formation (Klypin et al. 1999; Moore et al. 1999). 

Recent numerical simulations by Kravtsov et al. (2004) 
confirmed the significant overabundance of dwarf-like dark- 
matter haloes in the models in comparison to the relatively 
small number of observed dwarf galaxies around the Milky 
Way. They also included a simple model of star formation in 
dark-matter haloes which, for the first time, can reproduce both 
the circular velocity function as well as the radial distribu- 
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Fig. 15. Maps of the H i column density (left) and the heliocentric radial velocity (right) of M31 HVCs 1 and 2 with a higher angular resolution 
of about 30". The contours in the column density maps are drawn from 1 • lO'^ cm“^ in steps of 3 • lO'^ cm“^. The contours in the radial velocity 
maps are drawn in steps of 5 kms^'. Both clouds show compact sub-structure which is not resolved in the 2' maps. 


tion of dwarf galaxies around the Milky Way. For the central 
50 kpc around M31 the model of Kravtsov et al. (2004) pre¬ 
dicts only 2-5 dark-matter haloes with associated gas masses, 
Mg >10® Mo, that have not undergone significant internal star 
formation. They concur with the conclusion of Thilker et al. 
(2004) that only a subset of the HVCs found near M31 might 
represent primordial dark-matter haloes while the remaining 
clouds might be the result of tidal interaction. Our new imaging 
results support these conclusions. A recent complementary H i 
blind survey for HVCs around M31 with the Elfelsberg tele¬ 
scope (Westmeier et al. 2005a), however, yields no evidence 
for any further compact HVCs beyond the boundaries of the 
GBT map of Thilker et al. (2004), although extended, filamen¬ 
tary HI emission was found by Braun & Thilker (2004) even 


at larger projected distances from M31 with the WSRT. The 
Elfelsberg survey has a 3cr H i mass limit of about 6-10'^ Mq. It 
reaches out to a projected distance of about 130 kpc from M31, 
but it covers only a limited azimuthal range. The non-detection 
of additional compact H i clouds may be an indication that the 
thermodynamic considerations of Sternberg et al. (2002) may 
apply. Their calculations suggest that the H i mass associated 
with a particular dark mini-halo mass is a sensitive function of 
the ambient pressure of the circum-galactic environment. This 
would imply a substantial gradient of decreasing neutral gas 
mass with galacto-centric radius for a given dark halo mass. 
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Fig. 16. Maps of the H i column density (left) and the heliocentric radial velocity (right) of M3 1 HVCs 3 and 5 with a higher angular resolution 
of about 30". The contours in the column density maps are drawn from 1 • lO'^ cm“^ in steps of 3 • lO'^ cm“^. The contours in the radial velocity 
maps are drawn in steps of 5 kms“*. Both clouds break up into sub-clumps which are hardly resolved in the 2' maps. 


6. Summary and conclusions 

We have observed nine fields around M31 and one field near 
M33 in H i with the WSRT to study some of the HVCs discov¬ 
ered by Thilker et al. (2004) and Thilker et al. (in prep.) with 
arc-minute resolution. In two of the nine fields around M31 we 
could not detect the high-velocity gas, presumably because the 
emission is too diffuse and faint. In the remaining seven fields 
we identify about 16 individual HVCs with angular sizes of the 
order of 10', corresponding to absolute sizes of the order of 
1 kpc if a distance of 780 kpc (Stanek & Garnavich 1998) is 
assumed for M31. The observed H i masses are in the range of 
about a few times lO^* Mq to 6 ■ 10^ Mq. Under the assumption 
of gravitational stability only a very small fraction of the order 


of 1% of the total mass of the HVCs can be seen in neutral, 
atomic hydrogen, and additional mass components like ionised 
gas or dark matter have to be the dominant constituents. The 
application of the virial theorem is only of marginal suitabil¬ 
ity because it is subject to a number of restrictions which, in 
many cases, are not stringently fulfilled. The estimation of the 
dynamical masses of five of the HVCs, however, emphasises 
the necessity of a large fraction of undetected mass. 

Twelve of the HVCs are crowded in an area of only about 
1° X 1° at a projected separation of less than 15 kpc from the 
disk of M31. This remarkable complex of HVCs appears to 
consist of several filaments of clouds intersecting each other 
at various angles. Clouds along the apparent filaments have 
similar radial velocities. A detailed comparison suggests that 
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Table 2. Parameters for the determination of the dynamical masses 
of M31 HVCs 1, 2, 6, 8, and 10. Sv is the observed radial velocity 
gradient across the angular extent i? of each cloud as extracted from 
the one-dimensional cuts in Fig. 4, 6, and 9 and from the map in Fig. 5. 
denotes the lower limit for the derived dynamical mass. 


object 

(kms“^) 

& 

(arc sec) 

Miyn 

(10^ Mo) 

M31 HVC 1 

14 

400 

>86 

M31 HVC 2 

10 

400 

>44 

M31 HVC 6 

13 

100 

>19 

M31 HVC 8 

35 

540 

>730 

M31 HVC 10 

37 

540 

>810 


the HVC complex might consist of two distinct populations 
of clouds which can be distinguished by their different radial 
velocities, velocity gradients, and mean densities. Moreover, 
the HVC complex is partly overlapping with the giant stellar 
stream of M31 discovered by Ibata et al. (2001), suggesting a 
connection with the stellar stream and a tidal origin for some 
or most of these HVCs. The densest clumps within this region 
have large internal velocity gradients which are not generally 
associated with the orientation of nearby diffuse filaments, pos¬ 
sibly suggesting significant local concentrations of dark mass. 
Another HVC near M31 is located only half a degree south of 
NGC 205 and at an identical radial velocity, although the ev¬ 
idence for a physical association remains circumstantial. The 
remaining HVCs around M31 are compact and isolated without 
any obvious connection with M31 or one of its known satellite 
galaxies. 

The HVC near M33 is not isolated but appears to be con¬ 
nected with the H i disk of M33 by a faint gas bridge. Both the 
HVC and the bridge share a common radial velocity gradient. 
This configuration is again suggestive of a late stage tidal ori¬ 
gin, but the progenitor of the H i gas remains unknown. 

The modest number of isolated HVC detections seem to 
be in agreement with the numerical simulations of Kravtsov 
et al. (2004) which predict only 2-5 primordial dark-matter 
mini-haloes within the field of the GBT survey by Thilker et 
al. (2004). However, a complementary H i survey of M31 with 
the Effelsberg telescope by Westmeier et al. (2005a), reaching 
out to a projected distance of about 130 kpc from M31, reveals 
no further compact HVCs beyond the boundaries of the GBT 
map down to an Hi mass limit of about 6 ■ 10 "*Mq. The con¬ 
centration of HI detections at relatively small galacto-centric 
radii may be consistent with the calculations of Sternberg et al. 
(2002) which suggest a rapidly declining neutral fraction as¬ 
sociated with a given dark halo mass for a decreasing ambient 
pressure in a circum-galactic environment. 
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